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Abstract 

In this paper we present a proof for automatic 0(a) improvement in twisted mass lattice QCD 
at maximal twist, which uses only the symmetries of the leading part in the Symanzik effective 



f^ ■ action. In the process of the proof we clarify that the twist angle is dynamically determined by 

^^ , vacuum expectation values in the Symanzik theory. For maximal twist according to this definition, 

^D I we show that scaling violations of all quantities which have non-zero values in the continuum limit 

^ ' are even in a. In addition, using Wilson Chiral Perturbation Theory (WChPT), we investigate this 

Mh, definition for maximal twist and compare it to other definitions which were already employed in 



actual simulations. 



I. INTRODUCTION 

It becomes more and more apparent that twisted mass Lattice QCD (tmLQCD) [l|, |2| 
is a promising formulation to approach the chiral hmit of QCD, despite the fact that the 
flavor symmetry is explicitly broken. A twisted mass protects the Wilson-Dirac operator 
against small eigenvalues and therefore solves the problem of exceptional configurations, thus 
making numerical simulations with small quark masses feasible. Recent studies J3|, y , |5|, |6 1 in 
the quenched simulation were performed with mT^/mp values as small as 0.3 without running 
into problems due to exceptional configurations. Even though it will be challenging to reach 
such small pion masses in dynamical simulations, m.,^/mp < 0.5 seems fairly possible [7[. 
This numerical advantage of tniLQ_CD is supplemented by the property of automatic 0(a) 
improvement at maximal twist J8|j9|, |lO| . For a recent review of these and some more results 
in twisted mass LQCD see Ref. jll|. 

Some issues, however, remain to be fully understood. The proof of automatic 0(a) 
improvement in Ref. JS| makes use of the symmetries ruq = mo — nicr -^ —rriq and r — > — r, 
where ttiq is the bare (untwisted) quark mass, r is the parameter in the Wilson term and 
mcr(^) is the critical quark mass. Maximal twist is defined by setting the bare mass to a 
critical value, mo = m,cr{r). A concrete definition of m,cr{r) is not required in the proof 
as long as the symmetry property m,cr{—r) = —mcr{r) is satisfied, and tuning to the bare 
quark mass where the pion mass vanishes has been suggested as one particular choice for 
'^cr(^) (We call this definition "the pion mass definition" in the following). However, it 
has been pointed out in Ref. [ij] that the condition m.cr{—r) = — rracr(r) is violated for the 
pion mass definition by non-perturbative effects. Consequently, the O(a^) scaling violation, 
expected from an 0(a) improved theory, is lost unless the twisted quark mass satisfies the 
bound fi > a^Aq^j-,. Instead, terms linear in a and with fractional powers of a are predicted 
by Wilson Chiral Perturbation Theory (WChPT) for very small twisted quark masses. On 
the other hand, automatic 0(a) improvement is expected to hold if the critical mass is 
defined through the partially conserved axial vector Ward identity quark mass (PCAC mass 
definition). 

A .cent paper Q co.e. to a different eonc.n.on. Auto,nat,e 0«.) ..proven^ent Kas 
been proven by the m^ -^ —rrig symmetry only, without using the symmetry in r. It 
is claimed that both the pion and the PCAC mass definition guarantee automatic 0(a) 



improvement, but the remaining O(a^) effects differ significantly. In particular, the pion 
mass definition for mcr exhibits cut-off artifacts of 0{a^ /m?^) which are enhanced for small 
pion masses. These enhanced lattices artifacts are shown to be absent for the PCAC mass 
definition. 

12 1 and Ref. [l3[ are obviously in disagreement in the small quark 



The results in Ref. 



mass region where the bound /i > o^A^qy^ does not hold. Whether this is of relevance for 
present day simulations with lattice spacings a ~ 0.1 fm remains to be seen. A recent scaling 



analysis 



1J| of /tt in quenched tmLQCD seems consistent with an 0{a?) scaling violation, 



although its magnitude is rather large. Given the fact that automatic 0(a) improvement 
is a highly acclaimed feature of twisted mass LQCD, it is certainly desirable to study this 
issue further and find an explanation for these contradicting results. 



Close 



related is the so-called 'bending phenomenon', observed in quenched simulations 



|. The pion mass, the pion decay constant and the vector meson mass show an 
unexpected strong non-linear quark mass dependence for small quark masses if the pion mass 
definition for the critical quark mass is used. This curvature is significantly reduced when 
the untwisted quark mass is tuned according to the optimal choice proposed in |lj]. That 
this bending is indeed a lattice artifact of the twisted mass formulation is also supported by 
calculations using the overlap operator on the same gauge field configurations with similarly 
small pion masses. Here the bending is absent [3|. 

In this paper we revisit the property of automatic 0(a) improvement in twisted mass 
QCD. We first give an alternative proof for automatic 0(a) improvement at maximal twist 



without using the symmetries of t 
improved version of previous ones 



le parameters mq and r. Although our proof is just an 



a 



13|, we can clarify the meaning of "maximal twist" in 



the process of our proof. We will argue that, in presence of spontaneous symmetry breaking, 
the twist angle 9 is determined dynamically by the ratio of two vacuum expectation values 
in the Symanzik theory, namely 



te = ^ ' (V\ 

Provided that the mass parameters of the theory are tuned such that 6 = ±7r/2, we can show 
that the scaling violations of observables start with a^, i.e. the theory is 0(a) improved. 

We also investigate this new definition for maximal twist using Wilson Chiral Perturbation 
Theory (WChPT) 15|, |l6| (for a review see Ref. 17|). We explicitly show the absence of 



0(a, afi) contributions in the expressions for the pion mass and decay constant. 

We finally compare our new criterion with other definitions of maximal twist, the pion 
mass and the PCAC mass definition, which were previously employed in numerical simula- 
tions. We find that, although these two definitions show asymptotic a^ scaling violations, 
they do not exhibit the expected a^ scaling until a becomes small such that the bound 

^QCD 



/i > a^AgcQ is satisfied. 



II. ALTERNATIVE PROOF OF 0(a) IMPROVEMENT 
A. Main idea for 0(a) improvement 

The twisted mass lattice QCD action for the 2-fiavor theory is given by 

5'tmQCD = Sg + Stm- (2) 

The details of the gauge action Sq are irrelevant in the following, so we leave it unspecified. 
Stm denotes the 2-fiavor Wilson fermion action with a twisted mass term, which is defined 

as 

with 

[V+^i] (x) = ^ {U,{x)M^ + /i) - M^)) , (4) 

[v;^l] (x) = ^ (Mx) - ulix - fi)M^ - /.)) , (5) 

being the standard forward and backward difference operators. We supplemented the fields 
with the subscript " L" in order to highlight the fact that these fields are lattice fields. The 
parameters Mq and 6o denote the bare mass and bare twist angle. Instead of using this 
exponential notation it is also customary to write 

Moe''<^^'^' = mo + ^/io75r^ (6) 

where the bare untwisted mass tuq and the bare twisted mass /iq are given by 

mo = Mo cos 6*0, /io = Mosin6'o. (7) 

The lattice action (121) is invariant under the following global symmetry transformations 



1. U(l)(g> U(l) vector symmetry 

This transformation is part of the U(2) flavor symmetry of the untwisted theory. 

2. Extended parity symmetry 

PV : V^L ^ r^^^P^L, i'L ^ P^Lr'^\ (8) 

where P is the parity transformation, given by 

P^iix, t) = i'jo^L{-x, t), P^l{x, t) = -iijJL{-x, t)7o. 

For the gauge fields Pp is equal to the standard parity transformation. Note that 
ordinary parity P is not a symmetry unless it is combined with a flavor rotation in 
the 1 or 2 direction. Alternatively, one can also augment P with a sign change of the 
twisted mass term /xq, 

P = Px[^o^-H, (9) 

which is also a symmetry of the action. 

3. Charge conjugation symmetry 

C : i^L^x) -^ ijoj2i'L{x)'^ , ^l{x) -^ ~i^L{x)^ilol2, 
together with the charge conjugation transformation for the gauge fields, U{x, //) — »• 

u{x,fiy . 

Besides these symmetries the lattice action is also invariant under hypercubic lattice rota- 



tions and local gauge trans: 
According to Symanzik 



ormations. 



19|,|20|, the lattice theory can be described by an effective con- 
tinuum theory. The form of the effective action of this theory is restricted by locality and 
the symmetries of the underlying lattice theory. Taking into account the symmetries listed 
above one finds |l^ 

5eff = So + aS^ + a^S2 + ■■■ , (10) 



where the first two terms are given as 

So = 5'o,gaugc + [ d^x ^(x) U,D^ + M^e^^^^^'l V^(x), (11) 



Si = Ci f d^xi;{x)a^,F^,{x)^lj{x). (12) 

'S'o.gauge denotes the standard continuum gauge field action in terms of the gauge field tensor 
Fn^. The second term in 5*0 is the continuum twisted mass fermion action. The mass 
parameters are renormahzed masses, and we assume the renormalization scheme in (l[.^ It 
is worth mentioning that there is no "twisted" Pauh term if^'^^T^cT^vF^yil) present in Si, since 



2 



such a term violates the symmetry in eq. 0- 

In addition to the effective action we have to specify the direction of the chiral condensate, 
since chiral symmetry is spontaneously broken. From the fact that the direction of the chiral 
condensate is completely controlled by the direction of the symmetry breaking external field 
(i.e. the quark mass) in the continuum theory, we can take 

o L J /3a 

where \im.MR^o^^^"^v^oov{Mfj) ^ 0.^ Here the vacuum expectation value (VEV) is defined 
with respect to the continuum action Sq. To say it differently, the VEV ()13p defines the 
twist angle Q in the Symanzik theory. The above condensate is equivalent to 

(#)so = t;(Mfi)cose, (14) 

{i)i^^T^i))s^ = v{Mr) sin^. (15) 



^ Other choices for the renormahzed parameters are of course possible, but at the expense of additional 
terms in 5*1 of the effective action [13 . We also assume that use of the leading order equations of motion 
has been made in order to drop some 0(a) terms in Si. Without using the renormalization scheme in [ij 
and equations of motion there would be seven terms present in Si instead of only one |l8| . However, this 
larger number of terms would not alter the conclusion of this section. 

^ This property can also be derived from a different point of view. Since parity is conserved at /iq = in the 

lattice theory, ipj^r^a ■ Fip does not appear without /io- This argument can be extended to the case where 
the parity- flavor symmetry is spontaneously broken for a certain range of the untwisted mass Mq cos Oq 
|2l|,|22,|23,|2J,|23 in the lattice theory. In this case, the charged pions become massless Nambu-Goldstone 
bosons in the lattice theory, associated with this spontaneous symmetry breaking in the zero twisted mass 
limit. Therefore, it must also become massless in the Symanzik theory in the same limit. This fact also 
tells us that explicit parity-flavor breaking terms such as ipij5T^ cr f^,/ Fi^^^ip must be absent in the Symanzik 
theory without /io- 
^ The computation of this condensate follows standard arguments where one first considers the theory in a 
finite box with 4- volume V. See, for example, the appendix of Ref. |l5|. 

6 



We now want to argue that the choice 9 = 7t/2 (or — 7r/2) corresponds to maximal twist. 
In terms of the mass parameters this is equivalent to Mr = fin and mji = 0. In this case 
the action and the VEVs become 

Sq = 'S'cgaugc + d^x i! [tm^m + iMn'j^T^] iIj{x), (16) 

Si = Ci d^x ilj{x)af,^F^^{x)ip{x), (17) 

{^^)so = 0, (18) 

{iPil5r^4j)so = v{Mr), (19) 

It is easy to check that 5*0, the continuum part of the effective action is invariant under 

i/j -^ e™^^"''V, V ^ ^e™^5"''', (20) 

and therefore also under the Z2 subgroup Ti of this continuous transformation, defined by^ 

Ti^ = «75^V, Ti4) = i'i'j^T^. (21) 

Since T^ = 1 in the space of fermion number conserving operators, which contain equal 
numbers of ip and tp, the eigenvalues of Ti are 1 (Ti-even) or —1 (Ti-odd). The crucial 
observation is that the VEVs (fTHj) and (fT^ are also invariant under this transformation. 
The symmetry (j2(Jj) (and its discrete subgroup Ti) is not spontaneously broken, hence it is 
an exact symmetry of the continuum theory. The 0(a) term 

aSi = aCi / d'^x ijj{x)a^uFfj,u{x)ip{x), (22) 

on the other hand, is odd under Ti. Therefore non- vanishing physical observables, which 
must be even under Ti, can not have an 0(a) contribution, since the 0(a) term is odd under 
Ti and therefore must vanish identically. This is automatic 0(a) improvement at maximal 
twist. ^ Note that non-invariant, i.e. Ti-odd quantities, which vanish in the continuum limit, 
can have 0(a) contributions. 



^ A similar argument using this symmetry has been given independently by S. Sint [26j- 
^ This argument does not rely on our particular renormalization scheme and the use of the equations of 
motion. All possible terms in Si are Ti-odd once the continuum part is invariant under the transformation 



The above argument gives just the main idea of our proof for automatic 0(a) improve- 
ment, and we will give a detailed proof in the next subsection. However, one of the most 
important points of our analysis is that the condition for maximal twist and for automatic 
0(a) improvement is determined dynamically by the VEV (ipaipi) in the Symanzik theory. 
More explicitly, the symmetry (|2m) or (|2ip of the continuum theory must keep the VEV 
intact, so that the symmetry is not spontaneously broken. This condition seems natural, 
since the symmetry (f^ corresponds to a part of the exact vector symmetry in the con- 
tinuum QCD at maximal twist. Note that this symmetry refers to the vector symmetry in 
the so-called twisted basis J2, IS|- After rotating to the physical basis the theory is invariant 
under ordinary vector rotations in the 1 or 2 direction. However, for the proof of 0(a) 
improvement in the next subsection we prefer to stay in the twisted basis. 

B. Proof of 0(a) improvement 

Let us consider an arbitrary multi-local lattice operator O^^l {{x}), where {x} represents 
Xi,a:2, • ■ ■ ,Xn, d is the canonical dimension of the operator, t = 0, 1 and p = 0, 1 denote 
transformation properties under Ti and parity P: 

Ti : OZf{{x]) - {-iroZf{{x]\ (23) 

P- OZf{{SM) - (-imtt'({-^>0)- (24) 

Here we do not include the dimension coming from powers of the quark mass in the canonical 
dimension d of operators. For example, 

<t''(^) = J2^lIv;v-m^), (25) 

and in terms of these operators the lattice action is given by 

^tm = Y. [^lat''(^) - «^ <t''(^) + ^0 Ol^\x) + /io 0^^\X)] , (26) 

X 

with untwisted quark mass rriQ = Mq cos 6*0 and twisted quark mass fiQ = Mosin6'o. 



The lattice operator O^^l corresponds to a sum of continuum operators 0^"^"'" (^n: non- 
negative integer) in the Symanzik theory as 

oo 

oi:f ^ oir = J2 ^""-^ E E ^Ln..otf-\ (27) 

n=0 t„,p„ in 

where n is the canonical dimension of the continuum operator 0^"^"'" which consists of ip, 
ip, Afj_ and D^ only without any mass parameters, and 

Ti : Olf--{{x}) - (-l)*"0:r'"({a;}), (28) 

p. ol:^"^\{x,t}) ^ (-i)P"0*;;^"'"({-x,t}), (29) 

with tn,Pn = 0,1. Here we distinguish different operators with the same {tnPn,'^) by an 
index i^. To have a total dimension d in the expansion in Eq. ()27p . the coefficients c^^L „ j^ 
must be dimensionless: cfj^^^^^^ = cfJ^^^iJg^,\og{Aa),mqa, ^loa), where g^ is the bare 
gauge coupling constant, log(Aa) represents log-divergences of the lattice theory with some 
scale parameter A introduced in the Symanzik theory, and nig = rriQ — rricr is a subtracted 
quark mass with an additive mass counter term rricr, which will be specified later. Note 
that we consider possible power divergences of lattice operators by including operators with 
n = 0, 1, ■ ■ ■ ,d — 1 in the expansion. 

The following selection rules among these operators are crucial for our proof of automatic 
0(a) improvement: 

t + p + d = tn+Pn + n mod(2), (30) 

P + if 1^0 = Pn + (#/io)n mod(2), (31) 

where #yUo and (#/io)n denote the numbers of /io's in O^^l and Cj^p„,n,i„, respectively. The 
second equality can be easily proven by the invariance of the lattice action Q under the 
P = P X [fiQ ^ —fio] transformation, Eq. (jHl). To prove the first equality (jHOI), we introduce 
the following transformation: 



V', 



A..i-x)\ 

(32) 



A^{x) -^ -A^(-x)j 



which is a modified version of the transformation Vd introduced in Ref. J8[. Since it is easy 
to show that the lattice action (jS)) is invariant under Ti x P^, in addition to the invariance 
under Pp, the lattice action is invariant under Ti x V^ x Pp. On the other hand, combining 
the transformation property 

1 f (V+ + V,7) -/(x) ^ -(V+ + V-) -/(-x) 
\ [D,-f{x) - -D,-f{-x)) 

for an arbitrary function /(x) with eqs. (jH^ and (jH)), we can easily see that V^ x Pp counts 
the canonical dimension times the parity of the operator as 

^^ X P^ OZf{{x, t}) -. i-ir^'Ol^fm -t}) (34) 

V'.xP'p: Olf-'-i{x, t}) -. (^-i)n-^P^olf--i{x, -t}). (35) 

Therefore, the invariance of the action under Ti x V^ x Pp implies the first equality (j3(Jj) . 

Let us show how these selection rules are used to determine the structure of operators 
in the Symanzik theory. As an example we consider the operator O^^^ (x) in eq. (j25p . Since 
t + p + d = 4, the first selection rule gives t„ + (i„ + ra = mod(2), which leads to 

„01,3 
^^01,3 _ '-00,0 ^^00,0 01,3x.^l0,3 , 01,3^.^01,3 I ^ 01,3 x.^00,4 01,3 x.^00,4 01,3^^11,4 

^cS — ~Z3~ "T "-10,3*-^ "T "-01,3*-^ "T " "^00,4,^*^^ "•" " '^00,4,S*^B "'""'^11,4*-^ 

I 2 01,3 /o01,5 , 2 01,3 /oOl.S , 2 01,3 /ol0,5 , 2 01,3 /ol0,5 , fia\ 

where 

000,0 ^ ^ (^01,3 ^ ^fji^^r^^^ 010-3 = ^^, 

O'a'' = ^il^r^D^, Ol''' = ^t^.T^cr^^F^,^, Of'' = ^D^, Of' = ^a^^F,,^. (37) 
Applying the second selection rule that p + #/io = 1 = p„ + (#Aio)n mod(2), we obtain 

01,3/ \ -01,3/ 2 2\ 01,3/ \ -01,3/ 2 2\ 

Coo,o(/^oa) = /ioacooo(/^oa ), ^0,31/^00) = /ioaqo,3(/^oa ), 

01,3 / \ _ -01,3 1' 2 2n 01,3 / \ _ -01,3 1' 2 2\ /qg^ 

'^00,4,A(B) vA'-oO' j — A'-OO' Cqq^^^^^^^ (^/iga j, '^10,5,A(_B) vA'-oO' j — /iQCt Cj^q g^^^^^ (^/iQa J, (^OOj 

where only the /ioa dependence is explicitly written, and the other c^^^^ „ j^'s are even func- 
tions of /ioa. We then finally have 

/nOl'3 _ ^-2 -01,3x^,00,0 , „01,3xr,01,3 , ^2 -01,3 x/^00,4 2 -01,3 xr,00,4 

I 2 01,3 x/-)01,5 I 2 01,3 ^/-)01,5 
+ '^ '^01, S,/!*-^/! + ^ ^01,5,3^^3 

I -01,3^-010,3 I 01,3/^11,4 I 3 -01,3 //il0,5 , 3 -01,3 //il0,5 , /Qn^ 

10 



where all dimensionless functions are even in fiQa. It is important to observe that all opera- 
tors with t = contain only even powers of a, while those with t = 1 have only odd powers 
of a. 

Repeating the analysis given above for all operators which appear in the lattice action 
(|26p and introducing renormalized quantities (see appendix^ for more details), we obtain 



'S'tmQCD -^ ScS — 'S'o + niqSm + 2_^ l.*^ "'^2n + Ct '^ S2n~l\ ' 



(40) 



n=l 



where 



So = d^x I ipR (7^D^ + iiiR-i^T^) ^Ijr{x) - -FI^^r{x) 



rUgSm = mRS^ruR =mR d'^xipRipR(x) 



oO 



'-'2n-l 



d^X 



d'^x 



E^5 

i 

E2 



00,2n+4 ^R.i 



00,2n+4, 



X) 



10,2n+3 ^R,i 



■ o 



10,2n+3/ 



[X 



i 



01,2n+3 



^iR ■ o\ 



01,2n+3, 
1 

ll,2n+2 " t^R ' ^R,i 



R,i 



,ll,2n+2/ 



X) 



X] 



(41) 
(42) 

(43) 
(44) 



Renormalized parameters are introduced as^ 



/io = Z^{gR,\og{Aa),mRa,fx%a^)fiR, mg = Zj{gR,\og{Aa),mRa,n%a^)mR, 



A/2 



9 = Z(j (5fR, log( Aa),mRa,/i^a )gR 



(45) 



where gRj^UR, fiR are kept constant and finite as a ^ 0. We also define renormalized fields 



as 



ijjR = zl!'^{g]i,log{Aa),mRa,fiRa'^)ip, A/^^r = zl!'^{gl,log{Aa),mRa, fiRa'^)Af,. (46) 
A subscript R in O^'f"'" means that the operators are expressed in terms of renormalized 



2 ^2\ 



fields, and therefore ^i„p„,„ = Z*^p„„(^^, log(Aa), mRa,/i^a' 

Similarly, applying the selection rules to an arbitrary operator O^^l (again we give more 
details in appendix [X)) , we obtain 



ozf ^ of/ = y: 



l=-d 



E 



tp,d yi^[t+l]p,d+l 

'^[t+l]p,d+l,i ^R,i 



E~tp,d /?-\^ 

^[t+l]p,d+l-l,i ^R 



[t+l]p.d+l-l 



,(47) 



^ Note that the renorniahzation differs from the one usually employed in the Symanzik improvement pro- 
gram. 
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where [t + l] = t + l mod(2), p = 1— p, and coefficients c^^p^n^i ^^^ c^^'p^n.i ^^^ even functions of 
^jia. Note here that, even though we use the same notations as in eq. f|27j) . these coefficients 
are functions of gj^, log(Aa), mjia and fi^a"^, therefore the functional forms are different from 
the original ones. Formula (j47|) tells us that, if the lattice operator has t = 0, operators with 
t„ = in the Symanzik expansion appear with even powers of a while those with tn = 1 are 
associated with odd powers of a. This is reversed in the case that the lattice operator has 
t = 1: operators with t„ = are multiplied by odd powers of a in the Symanzik expansion 
and those with t„ = 1 by even powers of a. 

In order to obtain a finite result in the continuum limit, we have to remove possible power 
divergences in the expansion (j47|) by subtracting lower dimensional lattice operators from the 
original operator Oiat, in addition to subtractions of log(Aa) divergences including a possible 
mixing among operators whose canonical dimension is same as the original operator. We 
denote such a renormalized and subtracted operator as Oiatsub, which corresponds to 



c?(rti,sub ^ of^%^,^, = or + E «'^K.'^'' (48) 

1=1 

.n,[t+l\,d+l _ \^ tp,d ,^[t+«b,d+« ST- ;;tp,d .^lt+l]p.d+l-l /,QX 

^R;tpd — Z-^'-[t+l]p,d+l,i^R,i '^ t^R Z_^^[t+l\p,d+l~l,i^R,i A^^J 



where c? + / in the short-hand notation O^j.^J^ represents a canonical dimension of the 
operator and \t + /] labels the transformation property under Ti : 

T. : ^Kr^ - {-^r'0t:^t^\ (50) 

We conclude that, in terms of this general description in the Symanzik theory, the max- 
imal twist condition corresponds to the property that the continuum theory is invariant 
under the Ti transformation. This condition then entails mn = 0, which we call exact in- 
variance condition. However, it can be relaxed to mji = 0(a), which we call weak invariance 
condition. Imposing either of these we find 

1 /" I exact 

{Of')s,,^m,s„., = ^ \ V^^nV^l,RVA,,ne'^+--'-- Of' = (51) 

^ ^ I 0(a) weak 

for an arbitrary continuum operator O^' which is odd under Ti. (In the operator formalism, 
this condition expresses the fact that the vacuum |0) of 5*0 + m/jS*™^ is invariant under Ti. 
TilO) = or 0(a).) 
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Assuming the maximal twist condition is satisfied, i.e. mji = 0{a) at least, we now 
consider the following correlation function: 



\tp,d 



{orAsu^.i{^})) 



z, 



lat 



V^lV^lVU e^-Qc° Oitt R,sub({^}) 



(52) 



where Zj^t is the partition function defined by (1) = 1. In terms of the Symanzik effective 
theory, this correlation function corresponds to 



where we define 

{0)s = ^ j ViPuV^nVA^^ne^ O. 
For simplicity, we first consider the m^ = case. In this case we have 

{oo "^ oo 

n=l J n=0 

where we define a^S^^^ to be the sum of the a" terms in eq. ()55|). For example, the first few 
terms are given as 



(53) 



(54) 



(55) 



1\2 



Sio) ^ 1^ 5(1) ^ si 5(2) = 5° + ^ 
Under the Ti transformation, they behave as 

Ti : 5(") -^ (-1)"5("\ 
By expanding both action and operator, we have 



(56) 



(57) 



■\tp,d 



{0'^U{^}))s.. = ( 



Ot\{x}) + Y.a^0^^^t^\{x}) 



^a".S(")) 



So 



n=0 



= {Of\{x]))s, + Y. «'^"(^K/^'({^})^^"^)^o- (58) 

l,n=0,l+n^O 

Since terms with t + / + n = odd in the above expansion vanish from the maximal twist 
condition ()5ip. terms with t + 1 + n = 2s remain as 

oo 2s— i 

s=l 1=0 

oo 



St40f\{x}))s, + J2 «'^-*i^f "'({^}, 9l log(Aa), /i^; /i^ja^) 



s=l 



(59) 
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where we define 

Fr~\{x},gl\og{^a),^,n■,^^W) = i^(oK/+'({x})5(2-*-'))5o, (60) 

1=0 

which is an analytic function for small /i|.a^ (the last argument). This expression tells us 
that 

This proves automatic 0(a) improvement at maximal twist that scaling violations of Ti 
invariant quantities are even functions of a: a^""*^^ contributions are completely absent, while 
Ti non-invariant quantities have only contributions odd in a and vanish in the continuum 
limit. This completes our proof of automatic 0(a) improvement at maximal twist. (Here 
0(a") (n > 1) represents contributions of the form a""'"'^[log(Aa)]^ with s, /c = 0, 1, 2, ■ ■ ■ ). 

Notice that this proof for automatic 0(a) improvement is not restricted to on-shell quan- 
tities, and the equation of motion is not required at all for the proof. It is also noted that 
the proof does not require /ir = 0: Automatic 0(a) improvement is realized also for the 
massive theory. 

If rriR = 0{a) (weak invariance), the proof goes through with just a little modification 
(see appendix IX|). In the special case that m^j is odd in a {niji = af{a^)), we obtain eq. (j^ 
with a little modification in F^, while in more general cases with rriR = 0(a) the result 
becomes 

' 0{a), t= 1 

C. Ambiguity of the maximal twist condition in the lattice theory 

In this subsection we consider the maximal twist condition in the lattice theory and 
discuss the possible ambiguities of it. 

In the Symanzik theory, maximal twist is uniquely defined by the condition that an 
arbitrary Ti non- invariant operator O*"^^'"' has a vanishing expectation value, 

{0'P'')so = 0. (63) 
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Provided this condition is fulfilled, the expectation values of all Ti-odd operators vanish. 
Hence the particular choice for O^^''^ is irrelevant, and in that sense the condition ()63|) is 
unique. In the lattice theory, however, the maximal twist condition, which may be defined 

by 

{Ol^ksJ = 0, (64) 

depends on the choice of the operator C^i^'^sub' ^^'^ ^^ therefore not unique. In order to 
discuss this we make the Symanzik expansion of (j64p . which gives (see appendix ^ for 
unexplained notation and the derivation) 

= (^eff,R,sub) = o-H^if^n; a^ m% niRa, /i^a^) + mnH^ifiR, a^ m% niRa, ii\a^). (65) 

The solution 77^™'^^™'^' to eq. (J65|l . provided it is unique"^, is of the form 771™'^^™^'! = afi^a^), 
due to the symmetry of eq. (jF)K|l under the transformation mji -^ —m,f> and a -^ —a. 
Therefore, according to the analysis in the previous subsection, scaling violations in Ti- 
invariant quantities are even in a. 

A different choice for the lattice operator in ()64j) leads to a different solution jfj^aximai _ 
af{a^), so that the difference between the two definitions is of 0(a): Atti^'^^'™'^' = a{f{a^) — 
f{a?)}. Note that a solution 772™^^™^' in general depends on /i/j, inherited from the /xr 
dependence of H^. 

Let us consider some examples for maximal twist in the lattice theory. A simple one is 
given by 

((^^)lat,R,sub) = 0- (66) 

Unfortunately, this definition is not very useful in practice, since the subtraction of power 
divergences necessary for {ipip) prevents a reliable determination of this VEV in the lattice 
theory. Instead one may take Oj^x.y) = y4j^(x)P"(y) or (9iat(x,y) = 9^74^(x)P°(y) (a = 
1, 2), as was done in Refs. 3|, 14 Il4| : 

(A^(a;)P"(y)) = or {d,A';,{x)P^{y)) = 0, (67) 

where A°^ and P°' denote the axial vector current and pseudo scalar density, respectively. 
Yet another choice is J30| 

{Al{x)P%y)) = 0. (68) 



^ This seems plausible at small enough a and ^r, since the solution is unique in the Symanzik theory. 
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Depending on the choice for the axial vector current, either the local or the conserved one, 
the conditions (jUTj) lead to a different definition of maximal twist. However, the difference 
will be again of 0(a). 

We close this subsection with a final comment. Any maximal twist condition in the 
lattice theory determines a value for the bare mass niQ as a function of the bare twisted 
mass fiQ. It has been suggested to tune the bare mass to its critical value mo = rricr where 
the pion mass vanishes in the untwisted theory. However, this condition is not related to Ti 
invariance. For example, contributions from excited states violate eq. ()67|) even at ?7i^ = 0. 
Consequently, the pion mass definition does not correspond to maximal twist according to 
the Ti invariance condition. 

III. WCHPT ANALYSIS FOR 0(a) IMPROVEMENT IN TMQCD 

According to our analysis in the Symanzik effective theory, maximal twist is determined 
by requiring Ti invariance of expectation values. For example, for maximal twist we could 
impose 

{iP^) = 0, (^275^-3^) = v{Mr) ^ 0. (69) 

n this section we study this condition in Wilson Chiral Perturbation Theory (WChPT) 

nnnn 

15L llfiL l27L l28l |. and check explicitly whether 0(a) improvement is indeed realized. We 
also compare to some other definitions of maximal twist, which have already been used in 
numerical simulations. 

Automatic 0(a) improvement has been studied before in WChPT for various definitions 
of the twist angle and also for different power countings, which are determined by the relative 
size between the quark masses and the lattice spacing |l2, |22l, 120, uM ■ O^^ analysis follows 
closely the one in Ref. [l2|- We work mainly in the regime where m and fx are of O(a^) unless 
stated otherwise.* It is in this regime where the phase structure of the theor y is determined 
by the competition between the mass term and lattice spacing artifacts [18|,|32], and where 



the differences o 
contrast to Ref. 



t 



le various maximal twist definitions start to become relevant Il2|- In 



Q. 



12| | we work at higher order and include the terms of 0(?7ia, fia,a^) in our 



Since the parameters m and fi in WChPT are renornialized parameters we drop the subscript "i?" in this 
section. 
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analysis. These terms, which were also included in Refs. 3l|, |33| , provide a nontrivial check 
for automatic 0(a) improvement, since they are odd in the lattice spacing and, according 
to our Symanzik analysis, should not contribute to observables. 



A. Chiral Lagrangian and power counting 



In terms of the SU(2) matrix-valued field S, which transforms under chiral transforma- 
tions as S ^ LTiR^ the chiral effective Lagrangian reads 

£2 £2 £2 

+ (^4 + 1^5/2) (a^St^^E) (atS + Sta) 

- {We + Ws/2){m^i: + S%)(atS + S+a) 
-(iy^ + iy^/2)(ats + Sta)2 

- W^iia^E + E^df - W^2{a^a){a^J: + S^a) - Wx{d^m + rh^a). 



The terms through O(a^) have been previously constructed 



the last line, which were also included in Refs. 



31 



(70) 



m\M. Two terms of O(a^) in 



are easily derived with the spurion 



fields in Ref. J23- Note that the 0{ma) term proportional to Wx does not depend on S. This 
term is usually neglected since it does not contribute to the pseudo scalar masses and decay 
constant. Here, however, we will need it since it gives a contribution to the condensates. 
The coefficients /, B are familiar low energy coefficients of continuum chiral perturbation 
theory 



371, while all the W^s are additional low-energy parameters associated with the 

n n 

nonzero lattice spacing contributions [la,|35[- As usual, angled brackets denote traces over 
the flavor indices and the short- hand notation 



m = 2B{m + i/iTs) = 2Bme 



/iuj^Ti 



a = 2Wo a 



(71) 



is used |2S|- The mass parameters m and /i denote the renormalized untwisted and twisted 
mass J2|, which are defined according to^ 



m = Zmimo — rric 



jj, = Zf,jJo. 



(72) 



The renormalization constants Z^, ^/j are related to the renormalization const_ant Za of the axial vector, 
Za = Zm/Z^^ which follows from the vector and axial vector Ward identities 
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Even though the critical mass rricr includes the additive shift proportional to 1/a, it does 
not include certain contributions coming from the 0(a, a^) terms in the chiral Lagrangian 
[15[. For example, the third term in the first line of (|7Uj) gives rise to an 0(a) shift in the 
critical mass.^° 

Our power counting is based on the assumption that m '^ a? |l2l 
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where m stands 



for both the untwisted and the twisted mass and for p^ (proper powers of Aqcd are, as usual 
in this type of argument, understood). Since m and a are smaller than one we have the 
inequahties 

m ^ a > ma ^ a > m ^ ma ~ a . (73) 

According to this power counting the terms of 0(m, a^) in the chiral Lagrangian are of 
leading order (LO), while the 0(ma,a^) contributions are of next to leading order (NLO). 
Note that the size of the 0(a) term does not matter for the power counting, since it only 
contributes to the critical quark mass. 

B. Gap equation 

Starting from the chiral Lagrangian a gap equation for the ground state of the chiral 
effective theory can be derived. From the NLO expression of the chiral Lagrangian we find 
the potential 

V^ = ^25m'(ptS + T)P) + ^—2WQa{l + c^a^){T. + S^) - {^C2a\T. + St)^ 

+ ^£2am'(S + St)(pts + Etp) + |!c3a3(S + St)3 + ^^£!^cxa(Pt + p). (74) 
lb d4 4 

Here we introduced P = expicu^rs with tancuj^ = /z/m, and the following combinations of 
low energy parameters: 

C2 = -32 {2WI, + W'^)^, 62 = 32 {2W6 + W^ ' ^°^ 



(2Wof Wl Wn 

J I J 



cs = 64iy,i^-^, c~3 = 32W,2^, cx = 8Wx^. (75) 



These parameters are dimensionful and have [02] = 4, [02] = [£3] = 2, [c^] = 5 and [cx] = 1. 
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This term is often absorbed in the imtwisted mass, giving rise to the so-called shifted mass |15ll30l l31| 
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Since a twisted mass term breaks flavor symmetry we make the ansatz 

So = e'^^^ (76) 

for the ground state, and this ground state is determined by dV^/dcp = with 



+ fc2a{mt + /iVl - t^)t + fc^aH^ + 2Bfcxma. (77) 

Taking the derivative with respect to in ()77|) we obtain a gap equation for 

t = COS0, (78) 

which can be brought into the form 



Vl-t^ [x-t + 2/3„t + 7t2] = a[t- /3^(1 - 2t')] , (79) 

where we introduced the dimensionless parameters 

2S/i _ 2Bm + 2Woa(l + C3a'^) 

2c2a^ ' 2c2a^ 

C2am C2a Scsa^ . . 

In the following we will assume 

\f3^\ < 1, \f3m\ < 1, |7l < 1, (81) 

which can be justifled by a naive dimensional analysis when all dimensionfuU constants are 
assumed to be of O(Aqcd) together with the conditions aAqco < 1 and m/AqcD < 1- 

Note the sign convention for the coefficient C2- A positive sign corresponds to the scenario 
with spontaneous parity-flavor breaking |l5|, which guarantees the existence a massless pion 
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21|. A negative coefficient C2 results in a scenario with a first order phase transition |l8l.l3 
The details of the discussion of 0(a) improvement differ depending on the scenario for the 
phase diagram. In the rest of this section we are mainly interested in the scenario with 
C2 > 0, where spontaneous parity-flavor breaking causes some subtleties for automatic 0(a) 
improvement. These subtleties are absent for C2 < and we come back to this scenario at 
the end of this section. 



^^ See also Ref. [SlJ where it has been shown that the NLO terms in the chiral Lagrangian do not change 
the existence of two quahtatively different scenarios for the phase diagram. 
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C. Condition for 0(a) improvement in WChPT 

Taking derivatives of V^ with respect to m and /i, the two VEVs {'ipip) and {ipi'jBT^ip) are 
easily computed with the result 

(^^) = p^ = 2fB[{l+(3^t)t + cxa] 
am 



iijz^.Tsij) = j^ = V'B{1 + f3,t)VT^¥. (82) 

Therefore, the Ti invariance condition (j69|) corresponds to t = —cxa + 0{a^) in WChPT. If 
general scalar and pseudoscalar operators are employed for ipip and ipij^r^ipy these results 
are modified by 0(a) stemming from the effective operators in the Symanzik effective theory 
|4o|. This leads to 

(#) = -A^[{l + csat)t + csa], 

(iPi^.Tsij) = %^{l + cpat)Vl^l^. (83) 

Zjp 

Nevertheless, even in this case the condition (j69p leads to a similar result: t = —csa + 0{a^). 
We would find a similar result using the alternative condition (A^P^) = 0, which is 
equivalent to cot cjwt = with 

cotoJwT = -jyTp2J- (84) 

Note that Al is Ti-odd while V^ and P^ are Ti- even. Provided Noether currents are used 
for Aj^ and V^ one finds (see also appendix IBJ) 



(0|Aj=ok2) = fm^t{l + coat), (0|V^Uk2) = /m^Vl - t^ (1 + coat), (85) 

where we defined the coefficient 

Co = 16(21^4 + 1^5)^, (86) 

in analogy to the definitions in Eq. (f7^ . This leads to cot^wT = cot0, hence the condition 
{A'tP'^) = implies t = 0. The result differs if one uses general non- Noether currents. 
Additional contributions of 0(a) appear in the effective operators in the Symanzik expansion, 
which carry over to the chiral effective theory as well: 

(01^2 IvTa) = ^\t{l + CAat)-CAa\, (0|KLok2) = ^VT^^ [1 + cy at] . (87) 

Za L J '^ Zy 
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Note that here the currents on the left hand side are bare currents, as one can infer from 
the exphcit appearance of the renormahzation constants Za^ Zy. The way we have written 
the expectation values correspond to what can be directly measured in a lattice simulation 
without the knowledge of Z^, Zy- For cotc^wT "we find 

t(l + CAat) - CAa Zy 1 . 

coto^wT = 7=== X ^ X TT 7' (^^ 

a/1 - t2 Za 1 + CyOt 

and the maximal twist condition cotc^wT = gives t = c^a + O(a^), which has the same 
form as in the case of the VEVs. Note here that imposing a non- vanishing value for cot cjwt 
is sensitive to the ratio Zy JZa as well as to ca, ca and cy. As a final example we consider 
the condition (Ajf^P'^) = introduced in pj]. Since 

(0|PV') = ^[t-(l-2t2)cpa] (89) 

we again find t = cpa + 0{a^). 

To summarize, imposing Ti invariance we find the condition t = Xa + 0{a?) with some 
constant X. This constant depends on the specific choice for the operator in the matrix 
element. Nevertheless, all definitions guarantee automatic 0(a) improvement, as we want 
to show next. 

It is instructive to first consider the simpler condition t = (which is equivalent to a 
vacuum angle = ±7r/2). In this case the pseudo scalar mass and decay constant of charged 
pions, m^ and f-,, are given by (see also appendix lB|l 

m^ = , ^— , (90) 

U = fVr^[l + coat]. (91) 

These expressions are valid for arbitrary t, but for t = they turn into the results familiar 
from leading order continuum ChPT, 

ml = 2Bfx, U = /• (92) 

Apparently there are no 0(a, a^) corrections in these results. In addition, the O(a^) cor- 
rections are also absent, but this is not as surprising as one might first think. The charged 
pions are the Goldstone bosons associated with the spontaneous breaking of flavor and par- 
ity in the theory without a twisted mass term. Hence they must become massless when one 
enters the broken phase, i.e. when fj, goes to zero. With the same argument one would also 
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conclude that no terms of order 0(a, a^) terms are present. The same argument, however, 
does not apply to the 0(/ia) terms, and their absence is indeed a non-trivial demonstration 
of automatic 0(a) improvement once Ti invariance is imposed. 

It is now simple to show that we can relax t = to the weaker condition t = 0(a) without 
loosing automatic 0(a) improvement.^^ Suppose that t = Xa with some constant X. If we 
insert this into ()90|) and (j91|) we find (after expanding the denominator) 



m„ 



2Bfi(l + 



C2 , X 



Xa^) +0(/ia^), (93) 



U = f(l + [Co - X/2]Xa') + O(a^). (94) 

Again no 0{a, fia) corrections appear. This demonstrates, within WChPT and at least 
for the two observables we have chosen, that t = 0(a), which follows from imposing Ti 
invariance, is a sufficient condition for automatic 0(a) improvement. 

D. Other conditions for maximal twist and 0(a) improvement 

In the following we want to compare the condition of Ti invariance to some other condi- 
tions for maximal twist which are proposed in the literature. In particular we are interested 
in definitions where the untwisted mass rriQ is set to a particular value and kept fixed as 
one varies the twisted mass /iq. Such definitions obviously have a practical advantage for 
numerical simulations. Finding the /io dependent value mo such that a matrix element like 
the ones in ()66p or (jUTj) vanishes is computationally quite demanding, in particular in dy- 
namical simulations. One can save a substantial amount of computer time if one does not 
need to do this tuning for each twisted mass one wants to simulate, but rather stay at one 
fixed value of rriQ. However, such definitions do violate Ti invariance for most /iq values and 
it is therefore not obvious how this affects automatic 0(a) improvement. This is the issue 
we want to study in this section. 

The following two definitions keep the untwisted mass constant and both have already 
been employed in quenched numerical simulations pl y, S ly| : 

1. PCAC mass definition. For a given (bare) twisted mass fiQ the untwisted mass ?Tio(/Uo) 



^^ A similar argument that the theory is 0(a) improved for t — 0(a) has also been given independently by 
S. Sharpe Q. 
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is first tuned such that the PCAC quark mass, defined by 

vanishes. Then mo(0) = hm^jQ^o'^ol/^o) is used as the choice for rriQ of the PCAC 
mass definition for all /iq. Therefore, this definition is independent of /xo-^^ 

2. Pion mass definition. The bare untwisted quark mass is set to its critical quark mass 
where the pion mass vanishes at yUo = 0. In practice this value is usually obtained in 
the untwisted theory by performing an extrapolation oimi data to the massless point. 

In order to study these two definitions we have to translate ("match") them to the corre- 
sponding ones in WChPT: 

1. PCAC definition: The denominator in Eq. (j^J3|) is found to be given by 

(OIPV2) = ^il+P,t) (96) 

so that the PCAC condition reads 

Zpt{l + CAat)-CAa 

mpcAc = /i— ^===-- = 0. (97) 

^A vl — t^(l + Coat) 

This leads to t = cacl + O(a^) for any non-zero fi. Keeping, for simplicity, only the 
leading term and setting t = cacl into the gap equation, one finds 

X = CAa[l - 2P„i - ICAa] (98) 

in WChPT. In the following we will assume this x to be smaller than one. This is in 
accordance with our previously made assumption that all dimensionful coefficients are 
of order Aqcd and aAgco < 1- 

2. Pion definition: We need the expression for the pion mass in the untwisted theory. We 
cannot simply take the /i — »■ limit of Eq. ()90p . since t is equal to 1 in this limit and 



^^ Note that one can choose mo(Mo) at a fixed non- vanishing value fj,o as a (different) PCAC mass definition. 
Employing such a definition requires the determination of mo(A'o) at only one /xq value. Taking the 
/io ^ limit, on the other hand, requires the determination of 77io(/j,o) for various twisted mass values 
and a subsequent extrapolation. Hence the latter is numerically more demanding. 
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the whole expression is ill-defined. Instead, we first use the gap equation and rewrite 
the pion mass as 



ml 



2c2a 



2 



VT^ 



I - 1 + 2/3„ + 7t + 2f3^a 



1 + Coat 
Here the limit /i — > is well-defined and the condition ?72^ = reads 



(99) 



X = l-2Pm-7. (100) 

In order to check whether 0(a) improvement is realized we have to verify that t is at 
least of 0(a). To do so we have to solve the gap equation ()79|) with the x values in (jnH|l 
and ()1()()|1 . respectively. It is not necessary to solve the gap equation exactly, approximate 
solutions will be sufficient for our purposes. 

Let us assume t <^ 1, since we are interested in the small t case. In this case we can 
neglect the t^ terms in ()79|) and obtain the approximate solution 

t^^^^i±V- (101) 

a + l-2Pm ^ ^ 



For the PCAC mass condition we set x to the value in ()98|) and find 



Here we rewrote (I^Hl) as 



C2/i + 2ca (c2 + C2W0/B) a^ 
-"" 2Bii + 2{c2 + C2Wo/B)a^- ^ ^ 



am = -^a2 + 0(a^), (103) 

B 

and dropped all but the leading term proportional to a^. Taking into account that the 
denominator in eq. ()102j) is always of O(a^) or larger for C2 + C2Wq/ B > 0, together with 
our convention that B > and /i > 0, t is always of 0(a) and our assumption that t <^ 1 
is consistently satisfied. ^^ The solution in eq. p()2|l is of the form t = aX which we used at 
the end of the previous section (cf. eqs. (j93|) and (j94j) ) . with X representing the fraction in 
()102j) . However, here the value of X does depend on the relative size between a and /i. For 



^^ Recall that we here consider the case with C2 > 0. If C2 + C2W0/B < 0, the assumption t <C 1 could be 
violated at some value of fi. Therefore, in the latter case, we exclude such values of /i in the following 
consideration. 
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small a and fixed /i such that 2B^ ^ 2 (c2 + C2W0/B) a?, we can expand the denominator 
and find 

t ^ a||, (104) 

so X = cijIB and t = 0(a). Hence our discussion in the last section can be applied and we 
find O(a^) scaling violations in physical observables. On the other hand, for larger a such 
that /i = O(a^) we expect a modification of the simple linear a dependence of t, and this 
leads to distortions of the expected O(a^) scaling violations. In the extreme case of small 
fixed \i and large a such that 2B\i -C 2 (c2 + C2IV0/-B) a^ we find 

t ^ acA, (105) 

as expected from the definition according to the PCAC definition. Even though we recover 
a constant X, it is different from ()104|) . Since the sign of the low energy constants are a 
priori not known, it is even possible that the slope of t changes sign, depending on the size 
of a. This is of potential danger when one analyzes numerical data assuming a simple O(a^) 
scaling violation. The non-trivial a dependence of the r.h.s. in eq. ()102|1 is likely to obscure 
automatic 0(a) improvement in the region where [i is of 0(a^).^^ 

Let us now turn to the pion mass definition. Inserting (jlOOj) into the approximate solution 
fjlOip we obtain 

^ C2/i + 2(C2 + cWo/i?)a , . 

'- ''25/i + 2 (C2 + £21^0/5) a^' ^ ^ 

In order to derive this result we rewrote pOO|) as in eq. ()103|) and used, for simplicity, only 
the leading term proportional to c? . As before we find t ~ ac2/2B = 0(a) for small values 
of a, and, since the denominator is the same as for the PCAC definition we expect the 
modifications to become visible once the lattice spacing is such that /i is of O(a^). The 
details of the modification will be different because the numerator differs compared to the 
PCAC mass definition. 

However, the crucial difference between the PCAC and the pion mass definition is that 
the approximation pOlj) will eventually break down for the pion mass definition, since t goes 



^^ We emphasize that the non-trivial behaviour in t = aX{a'^) does not mean that there are terms hnear in 
a in physical observables. 
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to 1 for a vanishing fi. In that case the t^ terms can no longer be ignored. Interestingly, 
the approximate solution p06|) gives the correct value t = 1 at U = even though this 



approximation can not be justified. A more careful analysis finds [l2| 

1 



t ~ 1 - (5, 6 2 



fi{2B + d2a) Y^^ r^A/^^^/^ 



O 



,2 



(107) 



_{c2 + d2Wo/B)a\ 

Therefore, the condition t = 0(a) for automatic 0(a) improvement is satisfied only for small 
lattice spacings where 2Bfi ^ 2 (c2 + C2W0/B) a?. Even though this bound is asymptotically 
satisfied as a —i> 0, at a given non- vanishing lattice spacing the scaling violation becomes 
sizable for small twisted quark masses. In particular in the region 2i?/i <^ 2{c2 + C2Wq/ B)a?, 
where t — > 1 in the /i — > limit, automatic 0(a) improvement fails. 

This failure is seen when one inserts ()l()7j) into expression ()9H) for the decay constant, 
for example. Ignoring the small correction coming 5"^ and higher powers we find 

U = /v^(l + Coa-coa<5), (108) 

which obviously has a term linear in a. We emphasize that the reason for the presence of 
terms linear in a is the leading 1 in t = 1 — 5, and not the correction 5 with the peculiar 
dependence on fractional powers of a, even though the overall factor y25 complicates the 
whole a dependence. The leading constant term is the crucial difference to the PCAC mass 
definition, where t = 0(a), and this constant term spoils automatic 0(a) improvement for 
lattice spacings where 2Bfi <^ 2(c2 + C2Wi)/B)a^. 

We summarize the results in this section as follows. Although the PCAC mass and pion 
mass definitions lead to 0(a) improvement for small enough lattice spacings, the asymptotic 
O(a^) behaviour can only be seen at lattice spacings where the bound fi ^ a^AQ,-;^, is 
realized for a given /i. If this bound is not satisfied, the naively expected scaling violation 
is compromised, in particular for the pion mass definition, but also for the PCAC mass 
definition, even though to a much lesser extent. Note that the bound 2Bfi ^ 2(c2 + 
C2Wq/ B)a^ excludes automatic 0(a) improvement for the massless theory in the case of the 
pion mass definition, which is not unexpected since it is identical to the massless (untwisted) 
Wilson theory. 
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E. Automatic 0(a) improvement for the C2 < case 

Since the condition for automatic 0(a) improvement discussed in the previous section does 
not depend on the details of the lattice QCD dynamics, it seems applicable quite generally. 
However, there are circumstances when conditions like (A^P"^) = or cotcjwT = cannot 

nn 

be satisfied. This is the case when the 1st order phase transition scenario of Refs. |15l . I42| is 
realized. 

Let us consider this case in WChPT. For simplicity we work at LO only and set C2 = C3 = 



in the 



bl 



15 



owing argument. In this case, if C2 < 0, a first order phase transition appears at 
32j | and t is given by |l2J | 



t 



-VT^^, a^<l, X^O- , a = -^. (109) 



fJ'n 



0, «^ > 1, X ^ 



Although the condition t = can be realized for x = 0, the twisted mass /i must satisfy the 
bound /i^ > /imin.5 where 

,2\ 2 



2 / 2c2a^ 



Therefore, automatic 0(a) improvement can only be realized for /i^ > /i^;^ , in contrast to 
the parity conservation definition for the C2 > case, where no restriction on /i needs to 
be imposed. Note, however, that the same restriction on /i (at LO) is required for the pion 
mass definition in the C2 > case. 

This argument does not change qualitatively when one includes the NLO terms, as has 
been done in Ref . [^ . The phase transition line is no longer a straight line in the niQ — jjq 
parameter plane. If the term with C2 is included the maximal twist condition which gives 
t = becomes fi dependent and reads x = — /iC2a/(2c2a^). Nevertheless, the conclusion that 
one has to stay above the phase transition line in order to be able to satisfy the maximal 
twist condition remains unchanged. 

IV. CONCLUSION 

In this paper we gave an alternative proof for automatic 0(a) improvement in twisted 
mass lattice QCD at maximal twist. Whereas previous proofs J8|,|l3| used symmetries of the 
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bare lattice theory such as niq -^ —mg and r — > — r, we have used only symmetries of the 
leading part of the Symanzik effective theory in our proof. A more important observation, 
however, is that a precise definition for the twist angle, and therefore a condition for maximal 
twist, is determined dynamically by the ratio of two vacuum expectation values in the 
Symanzik theory: 

At 6 = ±7r/2, which is equivalent to Ti invariance of the vacuum in the continuum theory, 
scaling violations for all quantities are shown to be even powers in a, as long as they are 
invariant under the Ti transformation. Non-invariant quantities, on the other hand, vanish 
as odd powers in a. It is also shown that the ambiguity for the maximal twist condition in 
the lattice theory does not spoil automatic 0(a) improvement. 

We also studied the Ti invariance condition in WChPT. As expected, for the pseudo 
scalar mass and the decay constant we find automatic 0(a) improvement. 

We finally compared the Ti invariance condition to two other definitions for maximal 
twist, the PCAC mass and the pion mass definition. Both definitions have already been used 
in numerical simulations. These definitions have the practical advantage that the untwisted 
mass is tuned to a fixed value independent of the twisted mass /i. Although both definitions 
give asymptotic a^ scaling violations for /i ^ q'^Aqqy), we have shown that the expected a^ 
scaling can be obscured once this bound is violated. Hence naive continuum extrapolations 
can be deceiving and may lead to wrong results for these definitions of maximal twist. Here 
a WChPT analysis seems to be indispensable for a controlled continuum extrapolation. 
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APPENDIX A: SOME DETAILS FOR THE PROOF OF 0(a) IMPROVEMENT 

1. Derivation of Symanzik action 

We apply the Symanzik expansion to all operators which appear in the lattice action, as 
was done to O^^^ in the main text. After a little algebra we obtain the following expression 
for the effective action: 

oo 
5'tmQCD -^ ScS = 5*0 + rUgSm + 2_^ l.*^ "'^2n + «■ " 'S'2„_iJ , (Al) 



n=l 



where 



and 



So = j d'x [ZfOT{x) + ZpZ.ii 0''^\x) + ZgO'^^'] , (A2) 

Sm = j d^xZF-ZmO^^'\x), (A3) 

SL = [d'xJ2 ^00,2n+4 Or''"^'(^) + E C'01,2n+3 -" O'^'^'^^^x) , (A4) 

•^ i i 

'S'2n-1 = d X /_^CiQ2n+3 ^i ' (^) + ^_^Cll,2n+2 /^ ^i ' (^)- (^5) 

"^ i i 

In the definitions of S*^", the superscripts t„ = 0, 1 represent the transformation property 
under T^. 

Ti : Si" -^ i-lY'^Si". (A6) 

All coefficients which appear in the expressions above, such as Z^F^^^rn,G} and Cjp,^, are 
dimensionless functions of g"^, log(Aa), niga and fx^a'^. They are given in terms of c^^'p^ nii 
but their explicit forms are unimportant except for ?Ticr, which is given by 

ry ry _ 10,5 00,4 2 2-01,3 /' A V^ 

Zjp Zjm TncrCl — C;^g3 — C]^q3 — //qO C]^q3 \^' j 

where /xa c^q^ = c^q^- Using the selection rule (j^ it is easy to show that c^q^, c^q^^ and 
C]^q'3 are even functions of /ioct . 

Using renormalized fields ()46p and parameters ()45p. we finally obtain ()40p - ()44p in the 
main text. 
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2. Symanzik expansion of operators 

Using the selection rules egs-llHUj) and (jHH), we here determine the structure of the 
Symanzik expansion for the lattice operator, given by 



n=0 tn,pn i 

In the case with d = 2s, the selection rule (J3U|) gives 



(A8) 



<^T = EEK'"M^S2M^f'^' + ^— -^' 



21 I Jp,2s /y~,tp,2l 
"tp,2l,i^i 



1=0 i 



2{l-s)+l j tp,2s ^^tp,2l+l , tp,2s ^^tp,2l+l 
" 1 %,2l+l,i^i "^ Hp,2l+l,i^i 



(A9) 



where t = 1 — t and p = 1 — p. Furthermore, using the second selection rule (|HT|). we have 



1=0 i 



'tp,2l+l,i^i ' ^^R^tp,2l,^^ 



Similarly, for d = 2s + 1 we obtain 



(AlO) 



oo 



1=0 i 



tp,2l,l i r"'^ tp,2l — l,l I 



(All) 



Combining the results and rewriting the operators in terms of renormalized fields we finally 
obtain eq.(|T7j). 



3. Expressions for mq = 0(a) case 

In the case of mji = 0(a), the expansion of e^''^ becomes 

f oo 

e^^ff = e^« exp i m^Sm^ + J^ [a^^S^ + a^'^-'Sl^^,] 



n=l 



oo I. 

m 



«* E TT^™.""^'"'' 



k,n=0 



Under the Ti transformation, it is easy to see that 



T, : St. -. i-l)'St,, ^(") -> (-1)"^("). 



rriR 



(A12) 



(A13) 
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Expanding both the action and the operator, and using the fact that terms with t+l + k+n = 
odd in the above expansion vanish by the maximal twist condition (|3T|. we obtain 

((^Slsub({^}))5.«. = SMOf')so 

oo 2s— t 

+ E E a''-''-'m'kF'/-'^\{^}. gl log(Aa), /i^; rrma, fiW)iAU) 

s=l k=0 

where 

2s-t-k 

k\ 

(A15) 



Ff-'^\{x},gl^og{^a\^^n■,mna,^^W) = ^ {0^!^X'({^})%'S^''"^"~'^)so 



is an analytic function for small mjia and /X/jO . This expression tells us that 

(C'effR,sub({a;}))s,ft = < ^ ^ (A16) 

[ 0(a), t = 1 

for nin = 0(a). 

If we take rriji odd in a such that mji = af{a^), we have 

oo 

+ ^ a2^-*Ff -*({a;}, ^^, log(Aa), /x^; a^, /x^a^) 






(0|'")s„ + a'F! + a^Fi + ■ ■ ■ , t = 
^ aF] + a^Fl + ■ ■ ■ , t = 1 

where 



(A17) 



2s-t 



i^r *({^}> al log(Aa), /i^; a^ /x^) = Y.(f(^'^^'^d"'''i{^}, dl log(Aa), /i«; aV^ /i^^a^)- 



k=0 



4. Maximal twist condition on the lattice 

The maximal twist condition on the lattice leads to 

oo 2s— 1 



2s-k~l^k ip2s~l,k 
1 k 



oo oo 



k=0 s=k+l 

aH^{fj.R] a^, m\, niRa, ji^a?) + mRH^{fiR] a^, m%, iriRa, fi^a^), (A18) 
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where 

oo oo 

TT&t 2 2 2 2\ _ V^ 2k Sr^ 2{s-k~l) ■^^2s-l,2k+S / 2 i /a N 2 2\ 

H^{fXR;a ,mji,mRa,Hjia ) = 2_^mji 2^ a ^ >F^ (gf^, log(Aa), /xr; mRa,/i^a ) 

fe=0 s=k+l 

for (5 = 0, 1, and we keep the dependency on gj^ and log(Aa) imphcit in H^. 

APPENDIX B: WARD-TAKAHASHI ANGLE, PION MASS AND DECAY CON- 
STANT IN WCHPT 

In this appendix we provide some details about the calculation of cot ujwt and the pseudo 
scalar masses and decay constant. At leading order in our power counting scheme this has 
already be done in Ref. [l2|, and we refer also to this reference. 

Our first observable is the twist angle tuwT defined in eq. dHU). Instead of eq. (jHH) the 
twist angle can also be expressed as 

The extra derivative gives rise to an additional factor of the pion mass in both the numerator 
and the denominator, which consequently cancels. 

The Noether's currents appearing in the correlators on the right hand side of eq. (jBljl 
are given by (a = 1, 2) 

rpB 






coa 



4 



1 + ^(S + St) , V,\ = -^{r^i^d,^^ + St9,S)), (B2) 



4 
4 



1 + ^(S + St)] , Al^ = ^(r,(S9,St - St9,S)). (B3) 

The factor involving c^a stems from the wave function renormalization due to the 0{p^a) 



Aa _ Aa 
^11 — ^0,fi 



contribution in chiral Lagrangian, cf. (fflUl . In Ref. (l2] cottuwT was computed without the 
0(a/i, a^) contributions. Repeating the calculation including these terms we find 

cot Ci;wT = , = cot d> . (B4) 



This is the same result as in Ref. [12] • The functional dependence of cotci;wT is unchanged, 
and the 0(a/i, a^) corrections contribute only indirectly through the gap equation. 

Result eq. (IB4I) assumes that the currents in the correlators are Noether's ones that stem 

ri 

from the vector and axial vector Ward-Takahashi identities [1|]. This means that the point 
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split currents must be used in the lattice simulation. However, local currents are often 
used instead of there point split counter parts. This introduces additional contributions 
proportional to a. Taking into account the leading corrections of 0(a) only we obtain 



fy T /a, local T/a 



Z 4 a, local A a 



l + ^a(S + St)], (B5) 



ca 



^PB^ 



a(S + S^) +^i_g^9^(S-S^), (B6) 



4 



4 

where cv,a and ca are additional coefficients parameterizing the lattice artifacts stemming 
from the currents. Using these currents the expression ()B4|) changes to 



Zv 1 (1 + acAt)t - CAa 

COtWwT = ^— 7 n -^ , B7) 

ZaI + acyt VI - t^ 



which is the result in Eq. 

In order to calculate the pion masses we expand S around the vacuum configuration Eq 
defined in fl7(j|l . We parametrize the field S in terms of the pion fields according to^^ 

3 

S(x) = S|exp(5^m,(a;)r,//)s|. (B8) 

1=1 

Using this form in expression ()74|) for the potential energy we expand in powers of the field 
IT. The contribution quadratic in vr leads to the pion mass formulae 

m^^ = 2 Bm cos {(p — CO l) + 2WQa cos (p — 2c2a^ cos^ (p 

+3c3a^ cos^ (p + 2c2am' cos (p cos{(p — col), (B9) 

m^g = m^^ + 2c2a^ sin^ — 6030'^ sin^0 cos — 2c2am' sin 0sin(0 — cUi). (BIO) 

Expressing this in terms of m and /i and taking into account the gap equation we can express 
the mass for the charged pion alternatively as 

"^TT^ = ; ^^- (Bllj 

Note that the expression (jBllll is not singular for t = 1. From the gap equation one can 
infer that t = 1 can be a solution only if /i = 0, and in that case the result ()B11|) is not 
well defined. For t 7^ one can, using the gap equation, rewrite the pion mass formula as in 
which is well behaved for t = 1. 



^^ This parameterization differs slightly from the one used in Ref. [12|. This difference does not affect any 
of the results in this reference. 
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The decay constant is conveniently computed with the so-called indirect method [l|, |43[ , 
which is based on an exact PCVC relation and does not require the computation of any 
renormalization constants: 



U = ^(OlP^lvr^). (B12) 

In order to calculate the decay constant according to Eq. ()B12|) we need the matrix element 
of the pseudo scalar between the vacuum and the one pion state, where the pseudo scalar 
in the effective theory is defined by (r± = ^^4p-) 

P± = Z!^ri + ^(s + St)l(r±(S-St)) n^ = ^^l^^ (B13) 



4i 



4 



V2 



The matrix element in eq. ()B12|) is readily calculated at tree level with the result 

(0|P±|7r^) = fBil+P,t). (B14) 

With the expression for the charged pion mass we obtain 

U = f{l + coat)Vr^ (B15) 

for the decay constant. 
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